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In scalar-tensor theories the scalar fields generically couple nontrivially to gravity. We study the 
observable properties of inflationary models with non-minimally coupled inflaton and Dirac-Born- 
Infeld (DBI) kinetic term. Within the assumptions of the priors of our Monte-Carlo simulations 
we find these models can generate new interesting observable signatures. Our discussion focuses on 
string theory inspired phenomenological models of relativistic D-brane inflation. While successful 
string theory constructions of ultra-violet DBI brane inflation remain elusive, we show that in 
suitable regions of the parameter space it is possible to use cosmological observables to probe the 
non-minimial coupling. Fortunately, the most observationally promising range of parameters include 
models yielding intermediate levels of non-gaussianity in the range consistent with WMAP 5-year 
data, and to be constrained further by the Planck satellite. 



1. INTRODUCTION 

There is substantial evidence indicating that the early 
Universe underwent a brief period of rapid, accelerated 
expansion, or inflation [TH3]. In addition to solving con- 
ceptual and theoretical problems of the Standard Big 
Bang model (such as the monopole, flatness and hori- 
zon problems), inflation provides a quantum origin for 
the seeds of large scale structures in the Universe and 
produces theoretical predictions that are supported by 
detailed measurements of temperature anisotriopies in 
the cosmic microwave background radiation (CMB) [4]. 
But, despite its many phenomenological successes, the 
microphysical origin of inflation remains unclear; indeed, 
ever since its inception, it has been a challenge to derive 
successful inflationary models from fundamental particle 
physics [5]- Fortunately, early work on D-brane inflation 
in string theory [6Tlll|. has fruitfully led to promising 
realizations of the scenario 1 . 

Many inflationary models are driven, in part, by an 
inflaton field with non-canonical kinetic term, such as k- 
inflation |14j . Ghost inflation [15] and Dirac-Born-Infeld 
(DBI) inflation [16 . In these models the speed of sound 
can be significantly smaller than unity. Because cosmo- 
logical perturbations travel at speeds less than that of 
light in these models, distinct observational signals are 
typically generated, including detectable levels of non- 
Gaussianity, e.g. [PoT - TO?] . As we demonstrate in this 
paper, non-Gaussianity is an additional observable pro- 
viding new information that can be used to distinguish 
between certain inflationary models. 

It is the nature of quantum field theory in curved 
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spacetime that scalar fields generically couple non- 
minimally to gravity |20j . In this paper we study ef- 
fects of gravitational couplings in inflationary models 
with non-standard kinetic terms. In particular, we are 
interested in determining how such couplings alter cos- 
mological observables. Our goal is to determine if it 
is possible to probe such couplings using measurements 
of non-Gaussianity together with measurements of other 
standard physical quantities. We ultimately focus on the 
theoretically well-motivated setting of a non-minimally 
coupled inflaton field in a DBI inflationary model; how- 
ever, the methods we develop are applicable to more gen- 
eral settings as discussed below. 

An outline of the paper is as follows: In Sj2j we provide 
a brief introduction to gravitational couplings in theo- 
ries with non-standard kinetic terms and DBI models of 
D-brane inflation. In Sj3j we derive analytic expectations 
for the observables of canonical and non-canonical models 
with conformal coupling. In fj4j we discuss the prospects 
of observationally discriminating between pure DBI the- 
ories and DBI theory with a conformal coupling term. In 
J|5j we describe the priors and general methodology used 
in our Monte-Carlo simulations. In fj6j we present our in- 
terpretation and an analytic understanding of the results 
of the analysis. Finally, we provide a summary of sig- 
nificant findings and further conclusions in £j7j Technical 
calculations are relegated to a series of Appendices. 

2. MOTIVATION 

The most general actions we would like to study are of 
the form 

S = Jd+x^g (f(R, <p) + P(x, p)) , (l) 

where X = —^g^d^ipd^ip, so that for a canonically nor- 
malized scalar field, P(X,ip) = X — V(ip). F is a series 
whose leading contribution is the Einstein-Hilbert term 
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and further terms encompass possible gravitational cou- 
pling and correction terms 2 : 



F(R,<p) = *<p* 



*4i: 



R 

ip 2 



(2) 



Gravitational couplings of the form ^ are commonly 
generated in dynamical 4D gravity, and have been ex- 
plicitly computed for D-brane probes in curved back- 
grounds [HJ HI] > and in D-brane inflationary models [IT] . 
For a recent discussion of the effects of loop corrections 
on slow-roll inflation models in a power- counting formal- 
ism of effective field theory see [53] . 



2.1. DBI Inflation 

From a theoretical perspective, a particularly well- 
motivated instance of non-standard kinetic term occurs 
in D-brane inflation: the inflaton kinetic term is of the 
Dirac-Born- Infold (DBI) form [16] : 



p(x, <p) = -/- »vi - 2f(<p)x+r i (<p)-v(< P ) , (3) 

where the warp factor / is related to the harmonic func- 
tion of a warped compactification. The non-standard 
form of the kinetic Lagrangian imposes an effective speed 
limit on the inflaton field, analogous to the speed limit 
imposed on particle motion in Special Relativity. This 
behavior results in a new form of slow-roll inflation even 
when the potential is very steep [HJ [TS] . 

The most well studied D-brane inflation scenario con- 
sists of a brane anti-brane pair embedded in the warped 
throat region of a conformally Calabi-Yau flux compact- 
ification in Type IIB string theory. The warping can 
be used to obtain a large hierarchy in scales between the 
electroweak and Planck scales. If the throat is the warped 
deformed conifold |24j . it is well approximated by simple 
Anti-de-Sitter (AdS) space far from the tip. In the case 
of ultra-violet (UV) brane inflation, a mobile D3 brane 
falls from the UV end of the warped throat towards the 
infra-red (IR) end of the throat where an anti-D3 brane 
is located. The position of the D3 brane in the extra 
dimensions acts as an inflaton field as the brane moves 
toward the anti-brane under the influence of an attractive 
Coulomb force. Inflation ends via a tachyonic instability 
when the brane crashes into the anti-brane. The formal 
construction is now well known and we do not review the 
details of the model here. The curious reader is referred 
to the recent reviews jTSJ [13] . The most successful ver- 
sions of this scenario are slow-roll inflationary models. 



2 We work with the reduced Planck mass, written in terms of the 
four-dimensional Newton's constant Gjv as, Mpi = I/^SttGn ~ 
2.4 X 10 18 . 



We briefly comment on the difficulties involved in build- 
ing rigorous models of UV DBI inflation within the above 
setting. 

In the state-of-the-art constructions of D-brane infla- 
tion the inflaton field range in Planck units is bounded 
(23E6] by 



Aip 



Ah 



< 



pi 



(4) 



where N n » I is an integer representing the flux sta- 
bilizing the throat. The relatively few models capable 
of satisfying this bound and producing observable lev- 
els of non-Gaussianity typically predict a blue spectral 
index and a negligible tensor component [27] . Finally, 
DBI constructions must overcome difficult backreaction 
constraints ^6 ll28ti30] . 

In this paper we are interested in the general obser- 
vational properties catalyzed by combining the DBI ac- 
tion with gravitational couplings, and work within a phc- 
nomenological setting inspired by the afore mentioned 
brane inflation scenario. It is our modest goal to ascer- 
tain the possibility of observing gravitational couplings 
in inflationary models with non-standard kinetic terms 
and not to build realistic particle physics models. While 
all of the presented solutions conform to the latest obser- 
vational constraints imposed by the WMAP satellite [3], 
the majority do not obey the constraint on the field range 
Q. Interesting models of inflation in string theory capa- 
ble of avoiding this bound have recently been proposed 
in [3"lTl3"4"] . From the effective field theory standpoint, the 
DBI action is of particular interest even divorced from a 
particular string theory model as it represents the mini- 
mal realization of non-linear Lorcntz invariance. 



2.2. Equations of Motion and Gravitational 
Coupling 

Variation of the action ([!]) using Q with respect to the 
metric g^ v leads to the generalized Einstein equation: 



Fr R^ - -g^F + g^DF R - V m V„F,r = 
dptpdvtp - g^ v {d^ipdaip + f^if)) 

+ {r 1 (<p)-V(ip))g tw , (5) 

where Fr = dF/dR. The equation of motion for the field 

ip is: 

v^^+r 3 ^ -1 -!) 

—r^figrd^d^ - v-f' = q, (6) 

where / = 8/ dip, and we have defined 

1 
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y /l + f(tp)gi"'d ll ipd v <p 



(7) 
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In this paper we focus on a particular well- motivated and 
simple form for the gravitational coupling and truncate 
F(R, ip) at the second term in (El . This term, ~^Rip 2 , 



plays a significant role in brane inflation since the field 
associated to the brane position is a conformally coupled 
scalar [TTJ [16j EI] • The term is ubiquitous in effective 
field theories of quantum fields in curved spacetime |20| . 
and is renormalizable by power counting arguments. 



3. ANALYTIC ANALYSIS 

3.1. Canonical Models 

The effects of non-minimal couplings of the form (,R<p 2 
were first investigated in the context of canonical infla- 
tion, described by the action, 



S = 




2 2 ^ 



(8) 



in which the kinetic term of the field is canonically nor- 
malized, ip 2 /2. From Q we see that introducing the 
non-minimal coupling term generates an effective Planck 
mass of the form 



M, 



eff 



For field values larger than the critical value, 



M 



<y5 c = 



PI 



(9) 



(10) 



the effective Newton's constant can become negative. 3 
In a cosmological background with non-minimally cou- 
pled scalar (and positive £) the anisotropic shear diverges 
without bound as ip approaches <p c [35] [3BJ . We will re- 
gard such large field values as outside the range of validity 
of our effective theory. In all of the following analysis we 
ensure that tp < ip c . 

The observables of such theories are typically com- 
puted in the conformally-related Einstein frame, in which 
the inflaton is minimally coupled to gravity. The Einstein 
frame action is obtained by performing the well known 
Weyl rescaling of the metric, 



— o 2 



(11) 



where VL 2 = 1 - £ 



, , , , followed by a field redefinition, 
1 jm pi 

tp — y (7 (see Appendix [A| so that the action becomes 



3 For the present discussion wc will always take £ > 0. Subse- 
quently, we will focus on the conformal value 5 = 1/6 (cf. dis- 



cussion at the end of 



J d^x^g | 



M 2 ~ 



zgTd^dya - V(a) 



2 T 

' (12) 

where V(a) = V(tp)£l 4 (v?)- While this conformally re- 
lated theory is physically distinct from the original Eq. 
([I]), the cosmological observables calculated in this frame 
are equivalent [37] . Therefore, we conclude that any non- 
minimally coupled canonical single-field theory is obser- 
vationally equivalent to some minimally coupled canoni- 
cal theory (within the range of validity of the transforma- 
tion (11)). This makes it impossible to observationally 
determine the nature of the gravitational inflaton cou- 
pling in canonical theories. The result is in agreement 
with the analysis of [38] . 



3.2. Non-Canonical Models 

We next consider the case of a non-canonical inflaton. 
In what follows, we study DBI inflation [THJ [T5] as a pro- 
totype non-canonical theory. The standard minimally 
coupled DBI action is 



S = 




(13) 



1 



V(6) 



In brane inflation f(9) is the warp factor of the com- 
pactified geometry and V{9) is the potential that arises 
from the presence of anti-D3-branes and other perturba- 
tive and non-perturbative sources. The action for single 
field DBI inflation with non-minimal coupling is 



S = 




1 



M 2 



(14) 



f^g^D^d^ - 1 



As in the case of the canonical theory, we can remove 
the explicit nonminimal coupling via the transformation 
(111. This results in the action 



S = d x \/—g 



M 



pi 



R + P(X, , 



where 



The corresponding functional form for P is 

2 A2 
i-2 



(15) 



(16) 



i-e 



M 2 



P(cf>,X) = 6^X 



M 2 , 



-r i (0) 1 /i-2x/ i-e 



' M 2 



(17) 

+ r\4>)-v{4>). 
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Therefore, in terms of a homogeneous field 4>(t), the La- 
grangian (15 1 becomes: 



C 



2 



R + g(([>,4>) + 



(18) 



(i - sji-ftf^n^^-v^) 



where 



A2 j,2 



m|j ' 



(19) 



The presence of the function g((f>, (f>) makes it impossible 
to fully duplicate the observables derived from ( 14 ) using 



observables derived from a standard minimally coupled 



DBI theory (13). Because minimally coupled DBI infla- 



tion makes specific observational predictions [27J EH] > it 
is possible that this departure from a pure DBI theory 
might lead to distinct observable physics. As a first step 
towards exploring this possibility, we examine Eq. ( 18 ) 
under various field redefinitions. 



3.2.1. Non-Relativistic Limit 

We first redefine the scalar field so that its kinetic term 
reduces to the canonical form in the non-relativistic li mit . 
f((j>)<p 2 <C 1. In the non-relativistic limit the action (18) 
is reduced to 



Sn.r 



X V — 9 



PI 



R 



1-^(1-60^ V 
^-g^d^d^ - ~ 



Thus, we define a new field $ c by 



(20) 



dt> c = 

so that 

S nr = / d^x^J^-g 



(21) 



M'r 



M 2 ~ 1 - 

_E1 R _ -g^d^ c d v ^ c - V($ c ) 



(22) 

where V{<$> c ) = F(0)/fi 4 (0). Taking £ = 1/6, we can 
integrate (21) to obtain the relation, 



pi 



%/fjtanh | — — 



(23) 



Since the minimally coupled DBI theory ( 13 ) reduces to 



the same non-relativistic form as (22 ), we expect no phys- 



ical differences between the theories in this limit. 



3.2.2. Ultra- Relativistic Limit 

We next consider a different field redefinition which 
brings Eq. ( 18 ) into physical agreement with minimally 



coupled DBI inflation when the field is highly relativistic. 
In this regime, where the expression inside the square 
root becomes small, the warp factor plays a dominant 
role in determining the dynamics. This suggests a field 
redefinition </>—>$ that brings the expression inside the 
square root into the form present in minimally coupled 
DBI, 1 - $ 2 /($). Thus, we define a new field, 



Mi 



pi 



V M P] J 



(24) 



so that 



where 



C DBI = i 



C = C DBI +g{^,X), (25) 
- VW, (26) 



1-A/1-2/X 



and 



(27) 



(28) 



Here, $ c in the potential V($> c ) must be expressed in 
terms of $ by eliminating ip from definitions of $ and 
$ c , i.e. y($ c ($)). For £ = 1/6 we have the relation 



, ( $ c \ . ( * 
tanh — = = sm — = 

\V6MpJ \V6Mpi 



(29) 



In the relativistic limit, 2fX w 1, and since 
6£tan 2 (V^>/Af P1 ) « 1 for $ ~ ©(lCT 1 ) and $ ~ 
O(Mpi) (parameter magnitudes of interest in this pa- 
per), g is small and Eq. (25) reduces to the minimally 
coupled DBI Lagrangian ( 13 ). We elaborate on this issue 
in§! 



3.2.3. Intermediate Regime 

Since the above two field redefinitions result in a con- 
formal theory that is almost of the same form as a mini- 
mally coupled DBI theory in the non-relativistic and rel- 
ativistic regimes, respectively, we expect both theories to 
be observationally degenerate in these limits. However, it 
is unclear what the nature of the theory given by Eq. ( 18 ) 



is in the intermediate regime. Might there be a region of 
parameter space for which this theory gives different pre- 
dictions from a minimally coupled DBI theory? If DBI is 
the presumed inflationary theory, any observational dif- 
ferences between the two theories could yield direct in- 
formation about the presence of gravitational couplings, 
in contrast to the case of canonical models. In the next 
section, we pursue this question in detail. 
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4. OBSERVATION ALLY DISCRIMINATING 
NONMINIMAL THEORIES 

4.1. Non-Gaussianity 

Non-canonical inflation models have the defining char- 
acteristic that inflaton fluctuations generally propagate 
at speeds less than that of light. As a result, the fluctua- 
tions are non-Gaussian, at levels significantly larger than 
the negligible levels produced in single field slow-roll in- 
flation |40) . One non-Gaussian statistic is the three-point 
function of the co-moving curvature perturbation, £, or 
its fourier transform, the bispectrum [5TJ 22] : 



(C(k 1 )C(k 2 )c(k 3 )) 



(30) 



The bispectrum is characterized by its amplitude, 
shape, and scale dependence as a function of comov- 
ing wavenumbers kx, fc 2 , and fc 3 . For the leading-order 
contributions to the non-Gaussianities in general non- 
canonical theories we define the estimators [T5], 4 



Inl 



81 [ci 
35 / 1 
108 { cf 



- 1 - 2A 



where 



A = 



XP,xx + \X 2 P^xxx 
Px + 2XP XX 



(31) 
(32) 

(33) 



with P the fluid pressure, 
speed is defined 



The hydrodynamical sound 



P 



x 



P 



X 



2XP 



,xx 



(34) 



These two contributions correspond to different shapes 
in Fourier space [HI [43] . In the case of DBI inflation, 
fjq L is identically zero, and so an experimental detection 
of the shape of the non-Gaussianities offers the possi- 
bility of distinguishing DBI inflation from more general 
non-canonical scenarios, in particular, its nonminimally 
coupled counterpart. However, both shapes are maxi- 
mized for equilateral configurations (fci = fc 2 = k^) and 
will be difficult to distinguish in practice unless there is a 
detection of high significance [H] . We will therefore not 
consider using the shape of non-Gaussianities as a dis- 
criminator. Nevertheless, in [3S], it was suggested that 
since the speed of sound in the non-minimal DBI models 
differs from that of pure DBI, it may be possible to dis- 
tinguish between the two theories through cosmological 
observables. 



Our sign conventions for /jvl are opposite to those of the WMAP 
team. 



4.2. Observables 

We examine the full set of observables for each model, 
including the amplitude of non-Gaussianities as well as 
power spectra. We model the power spectrum as a power- 
law + running, 



Pn(k) = A s 



1— n s 4- h aln 



(35) 



and include a tensor contribution. We adopt a Monte 
Carlo approach [271 EH HZ] ~ we stochastically generate 
large numbers of inflation models of both types (min- 
imal and nonminimal) and calculate their observables, 
(n s , a, r, Jnl)- We then compare the distribution of pre- 
dictions of each theory in the observable parameter space, 
and identify regions that are populated by one model and 
not the other. One can then discuss the plausibility of 
distinguishing these models with future experiments. 

General non-canonical inflation leads to a power spec- 
trum of primordial curvature perturbations of amplitude 

na, 



Pn(k) 



1 



H 2 



8tt 2 M^ c s e 



(36) 



and a spectrum of gravitational waves of amplitude, 

2 H 2 



(37) 



Here, c s is as defined in Eq. ( 34 ) , and 



XP 



x 



M 2 X H 2 ■ 



(38) 



Because the density perturbations travel at speed c s , Eq. 
(361 is to be evaluated at sound horizon crossing, c s k — 



aPt, for each comoving wavenumber k. However, because 
the tensor perturbations propagate at the speed of light, 
Eq. (37) is evaluated when k = aH, with the result that 



scalar and tensor perturbations on the same comoving 
scale k exit the horizon at different times. Taking this 
difference into account can be important [48H5T] . For a 
given scale k, the tensor-to-scalar ratio 



Ph 
Pn' 



(39) 



can be evaluated at any time between this mode's scalar 
crossing and tensor crossing, since the quantities in it are 
assumed constant over this time period. For the solutions 
studied in this analysis, we have verified that the differ- 
ence in horizon crossings is not important [H t ~ H s ). In 
this case we may use the familiar 

r = 16ec s . (40) 

Additionally, we define the scalar spectral index, 

d\nP n (k) 



1 



dhxk 



-2e — 77 



(41) 



G 




FIG. 1: Monte Carlo results for 60,000 models each of DBI with conformal coupling (left) and minimal coupling (right). 
Coloring denotes the value of n s . Observables are quoted at k — 0.002 ftMpc -1 . 



where 



The equations of motion for the conformal models are 



e 

iff' 

Cs 

c s H' 



(42) 
(43) 



The running of the spectral index, a = gpt, can be ob- 



tained from the above via the relation, 



dink y/2e 
= Af P1 V 



d(f> 



1 



(44) 



In the conventions we are using, the action Q (without 
potential), is conformally invariant for the special value 
£ = g. The presence of a general potential or a non- 
standard kinetic term (such as in DBI) will break this 
symmetry. When it is necessary for us to consider an 
actual numerical value for the non-minimal coupling pa- 
rameter, we shall take it to be £ = |. We may also 
casually refer to the case of non-minimal coupling as the 
conformal case. 



5. MONTE-CARLO ANALYSIS 



5.1. Equations of Motion 



We obtain equations of motion for the minimal models 



from the Lagrangian (13) 



= r 1 (7-i) + ^, 
= -(p+p) 



(45) 



obtained from Eq. (25 I 



(46) 



r l {l-l) + 9,xX + V, 

-(P + P) 
= -(j + g)X. 

where 7 = (1 - 2JX)- 1 / 2 . 

5.2. Priors on V and / 



Care should be taken to appropriately determine the 
priors on the functions V(9), V($), and /(#), /($)• We 
must ensure that any observable differences found be- 
tween these models is truly a result of physics and not 
a poor choice of priors for these functions. Note that 
Eqs. ( |45j ) are equivalent to Eqs. (46) when g = 0, 
V(6) ~ V($), and f{6) ~ /($). This is because the 
Lagrangian Eq. (p6| coincides with that of minimally 



coupled DBI in the relativistic limit, when the warp fac- 
tor, /, is important. Therefore, to minimize physical dif- 
ferences, we should choose the same priors for f{0) and 
/($). However, Eq. (26) is not equivalent to the min- 



imally coupled Lagrangian in the non- relativistic limit, 
when V is important. Therefore, wc cannot impose the 
prior V{9) -H> V(Q). Rather, wc must appeal to the La- 
grangian (20), which does coincide with minimally cou- 



pled DBI in the non-relativistic model. This suggests the 
prior V{9) O V($ c ). 



5.3. Parameter Values and Methodology 

With the proper priors in place, we now define the 
functional forms of the potentials and the warp factors. 
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For the potentials we consider the Taylor expansion, 



v(9) = y: 



n=0 ' 



2n! IM, 



V"(* c ) 



V„ 



^2n! V^pi 

n— 



'PI 



4>, 



(47) 



(48) 



truncated at n = 4. Fourth order is about the extent 
to which we ever expect to reliably reconstruct the in- 
flaton potential, and hence this form is motivated by 
phenomenology. We assume that the warp factors are 
of the AdS form, f(8) = \/6 4 , /($) = A/$ 4 . Wc 
draw the potential parameters from the ranges |V n /Mo|j 
| V^/Vb| € [10~ 5 ,10 5 ], where V ,V Q is determined by 
normalizing the spectra to the best-fit WMAP5 value, 
A s = 3.2 x 10~ 9 [4]. The warp factor is determined 
by drawing the rescaled warping from the range A/Vo, 
A/Vb € [10~ 5 ,10 3 ]. We consider random initial condi- 
tions for Oi (and <I>i) as well as the initial field values, 
Qi (and <&i). While the minimally coupled field is 
unbounded, the conformally coupled field m ust satisfy 

$j < ^/|7rMpi, corresponding to the bound (10). The 

Monte Carlo analysis is carried out for the minimal mod- 
els by numerically solving Eqs ( 45 ) , and for the conformal 
models by solving Eqs. 



(46) 



Each solution is evolved forward in time until suffi- 
cient inflation is obtained, here chosen to be N = 55 e- 
foldings. The number of e-foldings is obtained by solving 
the equation 



dN = -Hdt 



(49) 



along with each system. Models which yield insufficient 
inflation are rejected, as are those that evolve to a meta- 
stable minimum and eternally inflate. For solutions that 
reach N = 55, by cutting off the evolution at this point 
we are supposing that inflation ends via tachyonic insta- 
bility in some auxiliary field. 

The correspondence between the comoving scale, fc, 
and the number of cfolds of inflation, TV, is different for 
the scalar and tensor modes because of the difference in 
horizon crossing times. It is determined by solving 



d(aH) 
dt 



= (1 - e)aH 2 



(50) 



with k = aH for tensors and kc s = aH for scalars. 
However, as mentioned, the difference in crossing times 
can be neglected, and so we evaluate all observables at 
aH = 0.002/iMpc" 1 . 



CL (with running and tensors) [3]: —0.1 < a < 0.03, 
0.9 <n s < 1.1, -253 < /^£ U < 151 5 , and < r < 0.3. 
We depict the fNL~ot plane and indicate the value of n s 
with the color bar, for the conformal model (left panel) 
and the minimal model (right panel) . There is a clear dif- 
ference in the clustering of points between the two mod- 
els. For the conformal model, we note the high density 
of points with red spectral index (n s < 1) with negative 
running (—0.1 < a < 0) and moderate degree of non- 
Gaussianitiy (1 < /nl ;$ 60). If there are observable 
differences between the theories, it is in this intermediate 
range that we expect the difference to be present, as the 
priors are chosen to ensure that the physics is the same in 
the non-relativistic (low Jnl) and the relativistic (high 
/nl) limits. The lowest limit in this range compatible 
with the sensitivity expected from Planck for equilateral 
type non-Gaussianity is f^ 1 > 10 [53]. We will demon- 
strate that the most significant differences in models can 
be seen in the approximate range /nl G (20,40). 



In order to demonstrate the accuracy of our numeri- 
cal treatment, we verify that the minimal and conformal 
models do indeed give the same physics in the DBI and 
slow roll limits. In the left panel of Fig. 2, we present 
models in the n s -r plane with a large degree of non- 
Gaussianity: 180 < < 200. Conformal models 
are colored blue and minimal models colored green. The 
two distributions of points overlap, indicating that the 
models are observationally degenerate in this range. In 
the opposite limit, f^ 1 < 1, we again see agreement, 
in the right panel of Fig. 2. However, when we focus 
attention on the intermediate region, 20 < ff?^ 1 < 40, 
in Fig. 3, we discover the surprising result that the two 
models yield distinct observational predictions. We find 
that there is an observable region that is populated by 
only one of the two classes of models; the conformal mod- 
els can accommodate a redder spectral index for a given 
r and /jvl ■ Additionally, the width of this strip is of the 
same order as the l-er error expected from Planck (gray 
error bar), suggesting that these models may indeed be 
resolved in practice. As we now show, the breaking of 
observable degeneracy is because of the difference in the 
field value at N = 55 between the minimal and conformal 
models. 



6. RESULTS 

We present the results of the Monte Carlo in Fig. 1 
for models that lie within approximate WMAP5 95% 



To begin, consider the following expressions for the 



A recent analysis finds -435 < /° q " < 125 [52] 



FIG. 2: Monte Carlo results for conformal models (blue) and minimal models (green) satisfying /jvi < 1 (left panel) and 
180 < /jvl < 200 (right panel). There is no detectable difference between the models in these ranges. 



DBI parameters and sound speed, 




FIG. 3: Monte Carlo results for conformal models (blue) and 
minimal models (green) satisfying 20 < /jvl < 40. The con- 
tours denote WMAP5+SDSS 1- and 2-a marginalized con- 
straints. The gray error bar is that expected from Planck at 
1-cr. 
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(51) 



where we have taken 7 w 7. It is now a simple matter to 
generalize the Lyth bound of [53] (see Appendix [C]) . 

Using the expression for the spectral index, Eq. (J4TJ) , 
and making various approximations valid for the regions 
of interest (see Appendix [B| we obtain 



1 



ry/3f NL 
4 



'2r 
X/Mm 



(52) 



where \ is a general scalar field denoting either (mini- 
mal field) or <I> (conformal field). This result indicates 
that the difference seen in Fig. 3 is a result of the dif- 
ference between #55 imax and $55 imaX i the maximum field 
values corresponding to iV = 55. While this is only a 
rough approximation, it provides a reasonably accurate 
understanding of Fig. 3. By choosing Jnl = 40 and 
taking the corresponding values of #55 iBlax an d ^ss.max 
(^.max > #55, max) in Eq. (52|, we obtain the curves 
approximating the cut-off of each model in the n s -r plane 
(cf. Fig. 4). In order to justify our choice of constant 
Jnl in this expression, in Fig. 5 we present the same 
conformal models as in Fig. 3, but with the value of Jml 
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0.2 




■ -J - 





0.01 - 



0.001 




0.0001 r 



le-05 



le "% 



FIG. 4: Same models as in Fig. 3, with theoretical curves ob- 
tained from Eq. ( |52[ ). The upper red curve approximates the 
cut-off in minimal models, and the lower curve for conformal 
models. 




FIG. 5: Same conformal models as in Fig. 3, with the value 
of /jvl color-coded. 



color coded. Note that models with the same /nl tend 
to form bands in the n s -r plane, increasing in magnitude 
towards small n s . These bands are of the same form as 
the theoretical curves found in Fig. 4. 

It remains to understand why $s5 jmax > #55,max m 
general. In the region under consideration, along the red 
curves in Fig. 4, we derive the approximate relation 



fV 



24 



(53) 



in Appendix [D] This expression is valid for both minimal 
and conformal models. We compare a minimal model 
with priors (f(0), V(9)) and a conformal model with 



FIG. 6: Close up of the negligible r portion of the parameter 
space. 



(/($), V($ c )). The approximate relation Eq. (531 im- 
plies that 

/($)n$ c ) - W)V(9). (54) 

Noting that V is an increasing function, while fV is a 
decreasing function, <i> c > $ implies that $ tends to be 
larger than 9. 

As an illustration, let us consider the case where 



f(9) = \9-\ V{9) =V + V 4 9\ 



In this case we have 



Y± 

V 



(55) 
(56) 



and thus $ > 9. Therefore, $ tends to be larger than 
9. Combining this with the result Eq. (52 1, we can say 
that, for /jvx = 40, the minimum value of n s for a given 
r tends to be smaller for conformal models than for min- 
imal models. 



6.1. Negligible r Regime 

Comparing Fig. 3 and Fig. 5, our attention is drawn 
to another region of potential interest near the Harrison- 
Zeldovich, n s = 1, r = point. Fig. 6 is a zoom of this 
region. To investigate this region of the n s -r plane we 
are required to improve on our estimates of Appendix B 
We provide a refined analysis to 0{^~ r ) in Appendix E 



6.1.1. N on- Minimally Coupled Models 

We begin with a discussion of the non-minimally cou- 
pled (blue) solutions in Fig. 6. We note that there ap- 
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pear to be no blue solutions with a red spectral tilt and 
negligible tensor contribution r < .06. This is easily un- 



derstood from equation (E3). The last term in (E3) is 
positive since s is negative. This implies that the non- 
minimal models have blue spectra (n s > 1) in the limit 
r — > 0. Hence, any detection of Jnl in the range [20,40] 
together with a measure of red spectral index and no 
detection of r < .06 would effectively rule out all non- 
minimal DBI solutions (within our prior selection). 



6.1.2. Minimally Coupled Models 

We now turn our attention to the minimally coupled 
DBI models (green solutions) . This brings us to a discus- 
sion of the large white region in Fig. 6 and the apparent 
absence of solutions with red spectral index and with 
low tensor signal r < .06. As is evident from our large 
/jvl analysis in Fig. 2, allowing models with significantly 
larger /nl{> 40) would gradually fill in this empty re- 
gion. The only successful red spectrum solutions are a 
small number of minimally coupled models located near 
the middle of the Figure. From (E5 ) with r ~ we find 



1 



2s 



3/ 



NL 



(57) 



where, in this region we have s < and \s\ = 0(1) and 
9/Mpi — 0(1). Hence, the reddest the solutions can be 
are those with Jnl = 20, giving n s ~ .97. The spectrum 
along the boundary of the populated region turns from 
blue to red at r = 7% , where 



1 (2\s\M/M n 



3/. 



NL 



1(T 



(58) 



Hence, any detection of f^L in the range [20, 40] together 
with a measure of red spectral index below n s ~ .97 
and no detection of r < .06 would effectively rule out 
all minimally coupled DBI solutions (within our prior 
selection). 

Solutions in the bottom right corner typically have 
X/Mpi — 0(1O~ 2 ), and therefore have a blue spectral 
tilt since we can no longer ignore the second term in Eqs. 
pb, p5b. 



• There is a strong degeneracy of observables in mod- 
els with very large and models with undetectable 
levels of non-Gaussiantiy. 

This result is not particularly surprising given our an- 
alytic analysis of [j3] However, for an observationally 
interesting intermediate level of non-Gaussianity with 
}nl € [20, 40] we find that it is possible, in princi- 
ple, to distinguish between minimally coupled and non- 
minimally coupled DBI models. 6 Assuming a detection 
of intermediate levels of non-Gaussianity /jvl G [20,40] 
and the value of the non-minimal coupling parameter 
£ = 1/6 we find: 

• For a given tensor-to-scalar ratio r and /jvx, non- 
minimal DBI models tend to have a redder scalar 
spectral index n s than minimal DBI models. In 
other words, for a given scalar spectral index n s and 
tensor-to-scalar ratio r, Jnl tends to be smaller 
for non-minimal DBI models than for minimal DBI 
models. 

Focusing our attention to the small tensor signal range 
we find: 

• Any detection of a red spectral index n s below .97 
combined with no detection of a tensor signal at 
Planck sensitivity effectively would rule out both 
non-minimal and minimal DBI inflation models. 

• Any detection of a red spectrum and no detection 
of a tensor signal at Planck sensitivity would rule 
out non-minimally coupled DBI models. 

Finally, the goal of this paper was to initiate the study 
of observable effects of gravitational couplings in theories 
with non-standard kinetic terms. While we have focused 
on the particular example of conformal coupling and DBI 
kinetic term, our methods are applicable to more general 
systems represented by the action ([I]). We leave further 
study of such models for future research. 



7. CONCLUSIONS 



In this paper we studied the effects of a non-minimally 
coupled inflaton to the DBI action. We assumed an AdS 
form for the warp factor and studied a wide range of po- 
tentials and parameter values using a Monte-Carlo analy- 
sis. Our aim was to determine the prospects for observa- 
tionally distinguishing non-minimal coupling from stan- 
dard minimally coupled DBI by adopting a phenomeno- 
logical approach rather than to conform with the latest 
model building technology in string theory. Our first 
finding was that: 



Although our priors are carefully chosen to minimize the differ- 
ences between the models, our conclusions are dependent upon 
these priors and the assumptions discussed in detail elsewhere in 
this paper. 
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Appendix A: Conformal Transformation 



Under the transformation 



_ q2 



(Al) 



with £1=1 — £t§2-, the Ricci scalar transforms as 



(A2) 



where we have omitted a surface term. From Eq. (Ill, 
the kinetic term transforms as 



- lV-99^d^d^ = -iv/q^TV^cW. (A3) 



The conformal transformation thus results in a modified 
kinetic term, 



(A4) 
(A5) 



The kinetic term can be made canonical by introducing 
a new field, a, satisfying 



due 



(A6) 



In terms of this new field, we obtain the Einstein frame 
action, 



M 2 ~ 1 

■'■\ -.'/ 1 - ^ v d^d v a - V(a) } , 



where V(a) = U^fT 4 ^) 



(A7) 



From the expression for the spectral index, Eq. (41 1, 
we obtain 

- (n s - 1) - 4e = (77 - 2e) + s 

3g,x(7 2 - 1)7 



2(7 + 5,x)(7 3 +.9,A") 



+ 



7 2 - 1 
9,x 



7 



3 7 3 - 7 + 2g. 



x 



#7 2(7 + .9,x)( 7 3 +9,x) 



J_ 
H Hf 

3g,xj 



2H 1 2 



x A 



Hf 



x B. 



(Bl) 



where 

A = 

B 



1-7" 



1 



(1 + 7 1 5 ! x)(l+7 3 9,x) ?>9,xl 1-7 



1 



1 



| 7 -a + | 7 -3 ffjX 



3g,x SxX _/_ 

27 ' ff,jr ' /'X ' (l + 7- 1 5,x)(l+7- 3 .9,x)' 

(B2) 



where \ is a general scalar field denoting either 6* (mini- 
mal field) or $ (conformal field). Because of the numer- 
ical facts 

7 = O(10), g iX = 0(1), ^ = 0(1), ^ = 0(1) 

(B3) 

we can safely set 



A~l, B ~ 1, 



7 
i? 7 2 



(B4) 



Thus, 



-(n s -l)-4e ~ \d,xc s s--^j = ^g iX c s s + ^, (B5) 

where we have used the ansatz / oc x~ 4 . As another nu- 
merical fact, we know that the first term in the final ex- 
pression is small (g tX = 0(1), c s = 0(0.1), s = 0(0.1)). 
In order to rewrite the second term in terms of observ- 
ables (and x/Mpi), we can use the relation 



2A 



2c 



{h X ) H 2 x 2 7 + 3,x 
Therefore, for x > and x < 0) 



Mm 

X 



2c s e 



(X/Mp,) 2 ' 
(B6) 



1 



-4e 



X/Mpi ■ 



From the expression 



16c s e, 



(B7) 



(B8) 
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we obtain 



n.„ - 1 



2r 



4c s x/Mpi ' 
Finally, by using the formula 

f - 1 

JNL 



we find 



n. - 1 



3c 2 ' 



'2r 



4 x/M P1 ' 
Appendix C: Generalized Lyth Bound 



(B9) 



(BIO) 



(Bll) 



(511: 



The generalized Lyth bound is derived from ( 40 ) and 

(CI) 



A/ 2 , 



Thus, we obtain 



i + i~ 3 g,x \h dt 



X (N = 55) - X (N = 0) 



Mpi r X55 
271 ixo 



i + 7 g,x 

(l + 7~ 1 3,x) 3 
3 



1/4 



Xo 



(02) 



where diV = —Hdt. 



Appendix D: Derivation of Relation (53) 



3(7 + 5,x) _ SM^ff 2 _ 1 



By using Eqs. ( |51| , (|46| and ([19 ), we find 

[7-l + /^]+5,x 

(Dl) 



2 7 2 



7 2 - 1 



h-l + fV]+g tX - 



Thus, 



7 2 -l 
2 7 2 



3(7 + S,x) 



7.x 



-7+1. 



From 



7 



1 / l + 7~ 1 ff,x 

c s y i + 7~ 3 g,x ' 



16c s 



(D2) 



(D3) 



we can rewrite Eq. (|D2|) as 

/V - I ■<•; 



1 + 7 ff,x 

1 + 7 _1 s,x 

1 / i + 7~ 1 g,x | 
c s y l + 7~ 3 #,x 



24 / (l + 7~ 1 g,x) 3 g,x 
r V 1 + 7~ 3 S,x 2 



This result applies to conformal models. On the other 
hand, for minimal models we have 



9x24 1 / 1 \ 24 /— — 

'1= I-ttt— 

V ZjnlJ r 



f V = (l-ci)--- + l=[l- 



(D5) 

This means that, for minimal models, the microscopic 
quantity fV is expressed in terms of observables r and 
Inl- 

Now, for /jvx = 40, the numerical facts 



c s ~ 7 - 1 ~0.1, g, x = 0(l) 
lead us to a rough approximation 



fV 



24 



(D6) 



(D7) 



for small r. 



Appendix E: Corrections of Order 0(7 1 ) 
1. Non-Minimally Coupled Models 



In order to improve accuracy of the formula ( 52 ) we 
can include corrections first-order in 7 -1 . The result is 



71., - 1 ~ — 



r\/3f NL 
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\/27 
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$/Mpi 


3g,x 
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4V3/wl 
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Here, we have used 
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2 r y ' X g tX /'$' 



' • - ~ 1 ( 1 + 2 7 ; -' / v 



2c, e 
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x (l - ^-^x) • 



(E2) 



By using / oc $ 4 and <?.x = tan 2 ($/v / 6Mpi), we obtain 



1 



4 

3ff,x 



1. 



(D4) 



4V37a^ 



$/Mpi 
(2$/ v / 6M P1 ) 
sin(2$/A/6Mpi) 



(E3) 
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2. Minimally Coupled Models 

The corresponding relation for minimially coupled 
models can be obtained by setting g x = 0. However, 
we have to be careful about the last term in A shown in 



(B2|. Then we obtain 



1) - 4e = 



20 



1 



J_ 



i?7 3 1 - 7 - 2 
This leads to the following approximate relation 



n , - 1 ~ - 



e/Mvi 3/. 



NL 



(E4) 



(E5) 



Interestingly enough, minimally coupled models have red 
spectra (n s < 1) in the limit r — > since s is negative. 
Here, we have used numerical values 20 < /jvl < 40, 
|s| = 0(1) and 0/Mpx = 0(1). 
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